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Abstract 

We analyze the quark spectral function above and below the critical temperature for deconfine- 
ment and at finite momentum in quenched lattice QCD. It is found that the temporal quark cor- 
relation function in the deconfined phase near the critical temperature is well reproduced by a 
two-pole ansatz for the spectral function. The bare quark mass and momentum dependences of 
the spectral function are analyzed with this ansatz. In the chiral limit we find that even near 
the critical temperature the quark spectral function has two collective modes corresponding to 
the normal and plasmino excitations in the high temperature (T) limit. The pole mass of these 
modes at zero momentum, which should be identified to be the thermal mass of the quark, is 
approximately proportional to T in a rather wide range of T in the deconfined phase. 



1. Introduction 

After the advent of Relativistic Heavy Ion Collider (RHIC), much attention has been paid for 
the properties of matter near and above the critical temperature of QCD phase transition, T c . To 
understand the structure of the matter in this region, it is desirable to identify the basic degrees 
of freedom of the system and their excitation properties. In the present study, we analyze the 
spectral properties of quarks, which are one of the fundamental degrees of freedom of QCD, 
at finite temperature above and below T c using quenched lattice QCD with clover improved 
fermions QjJ, l2|] . 

At asymptotically high temperatures, one can calculate the quark propagator using perturba- 
tive techniques. It is known that the collective excitations of quarks in this limit develop a mass 
gap (thermal mass) that is proportional to gT, where g and T denote the gauge coupling and tem- 
perature, respectively [3]. Moreover, in this limit the number of poles in the quark propagator is 
doubled; in addition to the normal modes, which reduce to poles in the free particle propagator, 
plasmino modes appear. 

In order to see the properties of quark spectral function p(a>, p) in this limit, let us first con- 
sider the spectral function at zero momentum. In this case, the Dirac structure of p(o», p) can be 
decomposed using the projection operators L± = (1 ± y°)/2 as 

p(co, 0) = (co)L +7 ° + p M («)L_ r °. (1) 

While p^(cl>) for free quarks are given by p^(a>) = 6(a> + m), in the high temperature limit one 
obtains p^(to) = [6(a> - mj) + 6(a> + mj)]l2, where mj - gT j V6 is the thermal mass. The two 
poles at a> = ±rriT correspond to the normal and plasmino modes. 
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The Dirac structure of p(a>, p) in the high temperature limit is also decomposed as 



P {co, p) = p p + (co, p)P + (p)y° + p p (w, p)P-(p)y°, (2) 

with the projection operators P±(p) = (1 +y°/»-y))/2 and p = \p\. The spectral functions p ? ± (to,p) 
in the high temperature limit read 

pl(co,p) = Zi(p)6(cJ + Ei(p)) + Z 2 (p)5{u> + E 2 (p)) + p mai .(±u), p), (3) 

where p CO nt.(w, p) represents the contribution of the continuum taking non-zero values in the 
space-like region, and poles at E\{p) > and E2(p) > corresponds to the normal and plasmino 
modes, respectively. The dispersion relation of the plasmino has a minimum at nonzero p Jst]. 



2. Quark spectral function at zero momentum 

In this section, we analyze the quark spectral function for zero momentum but with finite bare 
quark mass on the lattice. We use gauge field ensembles which have been generated and used 
previously by the Bielefeld group to study screening masses and spectral functions Jj,l21- The 
quark spectral function in this case is decomposed into p^(ci>) as in Eq. £[). In order to extract the 
spectral function p^(w) from the lattice correlation function, we assume a simple ansatz for the 
shape of p^(w) including few fitting parameters. We found that the two-pole ansatz forp^(w), 

p™{io)=Z l 6{co-E x )+Z 2 6{LO + E 2 ), (4) 

reproduces the lattice correlator quite well, where Z\%, and E\ 2 > are fitting parameters deter- 
mined from the correlated fit. The pole at a> = -E 2 in Eq. (O corresponds to the plasmino mode. 
The success of two-pole ansatz for the quark correlation functions suggests that the positions of 
the poles of the quark propagator would be near the real axis at to = E\ and -E 2 with small imag- 
inary parts. Provided that the positions of poles of the propagator are gauge independent, this 
also indicates that our results on the fitting parameters E\ and E 2 have small gauge dependence. 

In Fig.[T] we show the dependence of E\, E 2 and Z 2 j(Z\ + Z 2 ) on the bare quark mass, m p , 
for T/T c = 1.25, 1.5, and 3 obtained from calculations on lattices of size 64 3 x 16. The figure 
shows that the ratio Z%j{Z\ + Z£) becomes larger with decreasing m p and eventually reaches 0.5 
irrespective of T. The numerical result for each T shows that E\ = E% is satisfied within statistical 
errors there. These results mean that the quark propagator is chirally symmetric at this point |0]. 
Moreover, p^(o>) at this point has the same form as the spectral function in the high temperature 
limit. We therefore define the thermal mass of the quark on the lattice as m?- = (E\ + E 2 )/2 at 
this point. One finds that the ratio mj/T is insensitive to T in the range analyzed in this work, 
although it becomes slightly larger with decreasing T, which would be in accordance with the 
expected parametric form at high temperature, mj ~ gT. 

Figure Q] also shows that the relative strength of the plasmino pole, Z 2 j(Z\ + Z 2 ), decreases 
with increasing values of the bare mass, m p . The spectral function p^ico) thus will eventually be 
dominated by a single-pole. This result agrees with the perturbative result in Yukawa models that 
P^(cl>) approaches the spectral function of free quarks as the bare quarks mass becomes larger 
fii. 

For T < T c , we found that the lattice correlator is concave in the log-scale plot. This behavior 
indicates that the positivity condition of p^(w) is violated below T c , In fact, we have checked that 
the two-pole ansatz Eq. (O gives unacceptably large ^ 2 /dof below T c , It is also found that the 
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quark correlator does not approach the chirally symmetric one even in the chiral limit, in contrast 
to the previous result for T > T c . These results for T < T c seem consistent with a naive picture 
that quark excitations are confined and the chiral symmetry is spontaneously broken below T c . 

3. Quark spectral function at finite momentum 

Next let us analyze the quark spectral function at finite momentum on lattices with size 
64 3 x 16 for T/T c = 1.5 and 3. Throughout this section we consider the quark propagator in the 
chiral limit. The quark propagator in the chiral limit is decomposed into p ? ± {u),p) according to 
Eq. (12). Following the same approach used in the previous section, we adopt the two-pole ansatz 

p p + (o>,p) =Zi6(co-E 1 )+Z 1 6(co + E 2 ), (5) 

and determine four parameters from a correlated fit. The ^-functions at a> = E\ and -E2 corre- 
spond to the normal and plasmino modes, respectively. We found that^ 2 /dof with this ansatz 
is always smaller than 1.5 for all momenta analyzed in this study. This result means that the 
two-pole ansatz again reproduces the lattice correlation function well. 

In Fig. [2] we show the momentum dependence of the fitting parameters E\, E%, and Zil(Z\ + 
Z2) for T/T c = 1.5 and 3. The horizontal axis represents the momentum on the lattice p = 
(1/a) sin pa. The figure shows that for large momentum Zzl(Z\ + Z2) rapidly decreases and E\ 
approaches the light cone. The spectral function at large momentum therefore approaches that 
of a free quark, consistent with the perturbative result. One also finds that E2 is always smaller 
than Ei, in contrast to the results in the previous section. This behavior qualitatively agrees with 
the behavior of poles in the high T limit yfl. One also observes from Fig. [2] that £2 enters the 
space-like region at high momentum. 

3 



12 3 4 

p/T 

Figure 2: Dependences of the fitting parameters E\ and Ei and the ratio Zi/(Zi + Z2) on the lattice momentum p = 
(l/a)sin(/M) for T/T c = 1.5 and 3. 

An interesting property of the quark propagator in the high temperature limit Eq. © is that 
the dispersion relation of the plasmino has a minimum at finite momentum. In Fig. [2] one sees 
that the value of E% at lowest non-zero momentum on our lattice, p m ; n = 2nT(N T IN lT ) - 1.5T, 
is slightly larger than that at zero momentum, and the existence of such a minimum is suggested 
but not yet confirmed in the present analysis. 

In summary, we analyzed the dependence of the quark spectral function on temperature T, 
bare quark mass m, and momentum p in quenched lattice QCD with Landau gauge fixing. Above 
T c , we found that the two-pole approximations for the spectral functions in the projected chan- 
nels, p^(aj) and p ? ± (io,p), can well reproduce the lattice correlation functions. Although further 
studies on the volume dependence is needed, this result indicates that the excitations of quarks 
have small decay width even near T c . Below T c , on the other hand, the two-pole ansatz fails 
completely. 

The lattice simulations presented in this work have been carried out using the cluster com- 
puters ARMINIUS@Paderborn, BEN@ECT* and BAM@Bielefeld. 
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